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ALrsTRACT 

A  technique  was  developed  for  adapting  the  side-scan  sonar  to  obtain 
sonar  maps  of  the  undersurface  of  sea  ice.  The  equipment  is  portable  and 
air-mobile  operations  were  made  from  single-engine  light  aircraft.  Graphic 
records  displaying  the  range  and  relative  scattering  strengths  at  48  kHz 
were  assembled  into  sonar  maps  that  display  the  location  and  shape  of 
underlice  features.  Data  were  taken  at  five  sites  in  the  Arctic*  ,  (1)  2  km 
from  Pt.  Barro\%  Alaska,  (2)  a  hydrohut  near  Fletcher's  Ice  Island  (T-3) , 

(3)  2  km  from  the  edge  of  T-3,  (4)  175  km  from  Pt.  Barrow,  Alaska,  (5)  the 
AIDJEX  main  camp.  The  data  indicate  that  for  pack  ice  there  are  two  dis¬ 
tinct  types  of  backscattering:  very  higjh  level  backscattering  from  well 
defined  underice  ridges  and  very  low  level  backscattering  from  between  the 
ridges.  The  higher  scattering  at  the  ridges  is  probably  due  to  the  increase 
in  roughness  and  the  tilting  of  the  average  plane  of  the  scattering  surface. 
Comparison  of  the  sonar  map  and  the  aerial  photograph  shows  that  most  sur¬ 
face  features  have  subsurface  expressions  and  their  relationships  can  be 
complex. 


Ln.  t  c_. 

v^o  develop  the  underwater  technology  required  to  measure  the  location 
and  shape  of  features  on  the  bottom  of  sea  ice  and  to  compare  the  bottom 
with  features  on  the  surface,1  and  to  study  the  scattering  of  sound  from 

A 

the  water-ice  interface. 

<  f\ 

INTRODUCTION 

Bottom-scanning  sonar  has  been  used  to  measure  the  geographical 
location  of  features  on  the  sea  floor  since  World  War  II.1  The  morphology 
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anc  sound  scattering  are  related  as  follows:  large  rocks,  gravel  beds, 
etc.,  scatter  signals;  over  smooth  bottoms,  relatively  little  energy  is 
scattered  back.2  Somewhat  independently,  sonar  bottom  reverberation 
measurements  have  shown  quantitative  differences  in  the  scattering  over 
different  kinds  of  bottom.3  Sonar  reverberation  measurements  have  been 
made  over  bottoms  described  as  mud,  sand,  and  rocks  and  the  scattering 
strength  of  mud  is  usually  much  lower  than  sand  and  rock.  The  two 
approaches  need  to  be  put  together.  Theoretical  studies  of  the  high  fre¬ 
quency  round  scattered  at  rough  surfaces  have  shown  the  back  scattered  to 
be  dependent  upon  the  material  of  the  bottom,  its  shape,  and  to  a  lesser 

extent  the  frequency  of  the  signal  (for  roughness  much  greater  than  the 
acoustic  wave  length).^ 

Side-scanning  sonar  records  of  the  bottom  of  sea  ice  may  yield  infor¬ 
mation  about  the  morphology  of  sea  ice  and  scattering  of  sound  from  the 
sea- ice  bottom. 

METHOD 

Figure  1  shows  the  sonar  geometry  beneath  the  ice.  The  sonar  trans¬ 
ducer  is  on  a  rotating  mount  so  that  it  can  be  lowered  through  a  23  cm 
hole  that  is  drilled  with  a  gasoline-driven  ice  auger.  The  support  assembly 
consists  of  1.9  m  aluminum  sections  that  can  be  rotated  at  the  surface. 

The  transducer  can  be  lowered  to  any  depth  and  its  direction  can  be  con¬ 
trolled  to  less  than  1°.  The  sonar  is  a  modified  Kelvin  Hughes  Transit 
Sonar  that  transmits  a  1  msec  ping  at  48  kHz.  The  beam  is  fan  shaped  with 
a  beamwidth  of  1%°  in  the  horizontal  and  51°  in  the  vertical.  Underice 
features,  such  as  the  ridge  in  Figure  3,  scatter  sound  back  to  the  trans¬ 
ducer.  The  range  and  relative  backscattering  level  are  displayed  on  an 
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Figure  1.  Under ice  Sonar 
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m tensity-modulated  graphic  recorder.  The  recorder  has  a  dynamic  range 
afoot  20  db  from  white  to  black.  A  time-variable  gain  is  used  to  com¬ 
pensate  for  the  signal  level  decrease  as  range  increases.  Beyond  25  m, 
we  were  able  to  adjust  the  time-vat  fable  gain  so  that  the  record  was 
marked  at  the  same  intensity  for  comparable  scattering  features  out  to 
the  maximum  range.  Once  it  is  adjusted,  normally  the  gain  is  the  same  at 
the  same  range  for  subsequent  transmissions.  Maximum  ranges  of  275  m  and 
550  m  are  obtained  with  this  system.  A  picture  of  the  scattering  charactex 
of  the  ice  bottom  is  made  by  rotating  the  transducer  in  increments.  The 
data  are  displayed  by  assembling  display  increments  in  a  polar  mosaic. 

inderice  sonar  measurements  were  tnade  at  5  sites,  Figure  2.  The  first 
near  Pt.  Barrow,  Alaska,  was  used  for  instrumentation  tests.  The  others 
w?re  Fletcher  Ice  island  (T-3) ,  pack  ice  and  the  AiDJEX  main  camp. 

An  example  of  a  sonar  map,  taken  under  pack  ice  about  175  km  from 
barrow,  Alaska  is  shown  in  Figure  3.  The  transducer  was  at  depth  8.8  m 
-rd  the  ice  thickness  was  1.5  in.  For  easier  visualization  of  the  data,  we 
assembled  the  sonar  scattering  mosaic  using  3°  sectors  (Fig.  4).  The  data 
were  taken  at  different  gain  steps  to  increase  the  dynamic  range  of  the 
display  and  because  we  experienced  changes  in  the  overall  receiver  gain 
during  the  experiment.  Temperature-  or  velocity-depth  profiles  were  not 
made,  but  other  measurements  taken  under  arctic  sea  ice  for  this  season 
suggest  a  nearly  iso-velocity  sound  structure  for  our  working  depths.5,6 
Since  the  transducer  depth  below  the  underice  surface  is  much  less  than 
the  scan  ranges  of  the  sonar,  the  slant  ranges  shown  by  the  map  are  approx¬ 
imately  equal  to  the  true  ranges  of  the  underice  features.  The  sonar  map 
then  represents  a  plan  view  of  underice  features  that  scatter  sound. 
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The  sonar  map  shown  in  Figure  4  shows  the  edge  of  Ice  Island  T-3  and 
the  sea  ice  of  Colby  Bay  that  remains  attached  to  the  ice  island. 

UNDERICE  MORPHOLOGY 

The  relationship  between  surface  features  and  underice  features  is 
studied  by  comparing  the  sonar  teaps  and  aerial  photographs.  Figure  3  shows 
that  most  of  the  surface  features  have  subsurface  expressions.  An  under¬ 
ice  ridge  will  generally  scatter  sound  only  from  the  side  facing  the  trans¬ 
ducer  and  the  other  side  will  be  in  shadow.  The  data  are  superimposed  on 
igure  4.  At  this  site,  the  width  of  the  scattering  side  varies  from  about 
1  to  3  times  the  width  of  the  surface  expression.  The  prominent  subsurface 
ridge  (55  m  range  at  0°  and  the  ridge  trend  is  130°  clockwise)  is  generally 
deeper  than  the  transducer  and  has  blocked  sound  from  reaching  most  of  the 
underice  surface  beyond  it.  The  relationship  between  underice  features 

and  surface  features  can  be  complex.  For  example,  in  Figure  4,  the  underice 

■"  -  * 

expression  of  the  farther  surface  feature  appears  to  be  approximately  sym¬ 
metrical  with  its  surface  ridge,  whereas  the  underlce  expression  of  the  near 
feature  appears  to  be  located  only  on  one  side  of  the  surface  ridge. 

Sonar  data  were  taken  in  from  a  hydrohut  on  T-3.  A  sonar  map  is 
shown  on  Figure  5.  We  did  a  number  of  experiments  to  test  underice  sonar 
techniques  such  as  using  a  vertical  transducer  for  depth  contouring/  The 
interference  patterns  of  a  pair  ^ c  transducers  were  used  to  make  iso-angle 
contours  of  the  underice  surface,8  but  we  did  not  obtain  identifiable  zebra 
stripes..  In  the  smooth  areas  of  the  underice  surface,  the  scattering  levels 
were  too  low  to  observe  the  zebra’  stripes  and  in  the  rough  areas,  the  ridges 
were  too  rough  to  follow  the  zebra  stripes;  We  made  measurements  through 
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Figure  4.  Relationship  of  Surface  Features  and  nnderice  Scattering  Based 
on  the  Sonar  Map  and  Aerial  Photograph.  Surveyed  control  points  indicate 
that  the  aerial  photograph  Has  negligible  distortion  In  this  area.  An 
underice  ridge  will  scatter  sound  only  from  the  side  facing  the  transducer, 
and  the  other  side  will  be  in  shadow.  The  underice  ridge  at  the  top 
appears  to  be  more  or  less  symmetrical  with  its  surface  express  ion  j  whereas, 
the  lower  under ice  ridge  appears  to  be  a  projection  located  only  on  one 
side  of  its  surface  expression. 
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Figure  5.  Sonar  map  showing  edge  of  Ice  island  (T-3)  and  the  sea  ice  of 
Colby  Bay  that  remains  attached  to  the  ice  island.  The  edge  of  T-3  is  the 
strong  reflection  on  the  right  side  of  the  map. 
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the  hole  in  a  hydrohut  near  T-3  with  a  crossed  pair  of  transducers  to  ob¬ 
tain  a  pencil  beam. 

Most  of  our  depths  were  measured  by  placing  the  transducer  in  the  ver¬ 
tical  position.  We  found  that  this  technique  has  problems  with  resolution 
because  of  the  wide  horizontal  beam  and  the  vertical  side  lobes.  The 
results  are  useful  but  we  learned  that  it  will  be  desirable  to  shade  the 
transducers  to  reduce  the  side,  lobes.  Additional  information  was  ob  ained 
by  studying  the  changes  of  the  reflections  and  shadows  on  sonar  maps  taken 
at  different  depths. 

Both  sonar  and  surface  surveys  were  carried  out  at  site  #3.  Seven 
sets  of  data  were  taken.  Three  sets  have  a  narrow  horizontal  beam  and  four 
of  them  have  a  narrow  vertical  beaut.  On  Figure  6,  a  surface  map  superimposed 
on  a  sonar  map  is  shown.  We  have  estimated  the  subsurface  elevations  of 
the  features  by  using  the  information  collected  from  those  four  sets  of 
narrow-vertical -beam  data. 

The  surface  survey  included  a  survey  across  a  ridge.  We  have  shown 
on  Figure  7  the  under ice  scattering  features,  the  survey  path  and  the 
corresponding  cross-section.  The  estimated  sail  heights,  keel  depths  and 
the  ratio  keel  depth/sail  heights  for  the  four  ridges  shorn  on  Figure  7 
are  listed  in  Table  X.  The  average  of  the  ratio  is  7.1. 

Tab  In  t 


Estimated  Pen'*. 
Too  Elevation 

Estimated 
Keel  Depth 

Ratio  -feel  dgEth 
hex eh t 

1.425  m 

10  m 

7 

1.17  m 

9  m 

- 1 — 

7.69 

1.5  m 

11  m 

7.25 

1.48  m 

10  m 

6.75 

: 

*  &'■  ■ 

**  y-  v.  ■  .{IV 

3k  -  ,;rir 


,  ^v'-:  s 


lifc.  ■-  '■  1 1  |  v 

.  :<«r  ^ 


\4  M  vf!  •'•  /  v 


*  Y®  •  !:• 

}  i  r*  j> 


♦ 

% 


v 


■{f  i1/,  i* 

i  J  4 


It! 


»•*  I  : 

l  '•  \  •  I 


I  • 

*  n 
V;  m 

I- 


.  .v 
>i'£'S' ' 


& 


J;  ,.•--- 

ms#*-* 

rid!#*#*  - 


•■’;*  vV- v 


V 


*  / 


v  .,  ;■  .  — 

n. 


••■t* 

vi-' 


& 


SURFACE 

FEATURES 


500  i 


Figure  6.  Sonar  map  and  surface  features  on  pack  ice  near  Fletcher's  Ice 
Island  (T-3),  April  1972. 


Figure  7.  Scattering  features  and  profile  of  pack  ice.  The  surface 
!leVHt*on8  wef®  measured  along  the  survey  path.  Underice  depth3  are  from 
Tro  dotted  portions  are  in  shadowed  regions  for  the  sonar. 


sonar  data. 
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The  University  of  Washington  Applied  Physics  Laboratory's  subsurface 
vessel  (UARS)  surveyed  the  ice  under  Colby  Bay  area,  and  part  of  its  four 
tracks  extended  into  our  survey  area.  Their  results  agree  well  with  ours 
(Robert  Francois,  written  communication,  1972). 

The  morphology  of  the  undersurface  of  the  sea  ice  of  site  #4,  and  the 
relationship  with  the  top  surface  are  discussed  by  Berkson,  Clay  and  Kan, 
1973,  Journal  of  the  Acoustical  Society  of  America,  in  press.  Sites  #2 
and  #3  were  discussed  by  Berkson,  Clay  and  Kan,  1972,  American  Geophysical 
Union  Fall  Meeting.  Sites  #3  and  #5  are  discussed  by  Kan  (M.S.  thesis, 
1973);  Kan,  Clay  and  Berkson  (in  preparation,  1973);  and  Clay,  Kan  and 
Berkson  (in  preparation,  1973). 


SCATTERING  OF  SOUND 

t  -i .  - 

The  data  show  very  low  backs cattering  levels  at  48  kHz  of  the  under¬ 
ice  surface  between  the  ridges  and  very  high  levels  at  the  ridges.  The 
relatively  flat  underice  areas  between  the  transducer  and  the  first  ridge 
will  always  be  illuminated;  but  beyond  the  first  underice  ridge,  areas  may 
be  in  shadow.  The  length  of  the  shadow  depends  on  the  downward  projection 
and  the  transducer  distance  and  depth.  In  Figure  3,  the  relatively  flat 
areas  between  the  transducer  and  the  first  ridges  are  characterized  by  very 


low  backscattering.  (The  scattering  near  the  transducer  may  be  large  but 
the  rapid  change  of  the  time-variable  gain  precludes  qualitative  estimates 
of  scattering  for  ranges  less  than  25  m.)  The  very  low  values  are  indi¬ 
cated  by  the  lack  of  any  contrast  between  the  backscattering  background  of 
the  underice  surface  between  the  ridges  and  the  shadow  zones  behind  promi¬ 


nent  unde rice  ridges,  even  at  the  highest  |aln  steps.  However,  there  is 
very  high  backscattering  from  well-defined  areas  of  the  underice  surface 
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which  generally  have  surface  expressions.  Qualitatively,  the  scattering 
cross  sections  at  ridges  are  at  least  20  db  greater  than  the  flat  areas. 

In  the  conventional  backscattering  experiment,  the  backscatterlng 
levels  vs.  time  or  angle  are  measured  and  the  backscattering  function 
(scattering  cross  section  vs.  angle  of  incidence  or  grazing  angle)  is  de¬ 
rived  from  the  gebmetry  of  the  measurement  and  ray  tracing.  If  an  omni¬ 
directional  source  and  receiver  are  used,  one  usually  assumes  that  the 
scattering  interface  has  statistical  spatial  stationarity  for  data  analysis.^ 
For  example,  if  the  scattering  is  from  features  of  the  order  of  a  few  cm, 
then  any  area  of  the  order  of  a  few  square  meters  should  have  the  same 
roughness  as  any  other  area.  Then,  the  observed  backscattering  function  may 
be  compared  to  theoretical  functions  for  rough  surfaces.  Let  us  consider 
another  model  consisting  of  smooth  surfaces  surrounding  rough  patches  at 
random  distances.  In  omnidirectional  backscattering  experiments,  the 
nature  of  che  scattering  surface  would  not  be  revealed  because  the  addition 

i  ■ 

of  the  backscattered  signals  coming  from  all  directions  could  give  rela¬ 
tively  smooth  backscattering  functions  such  as  those  observed  in  measure- 

5  9  10  11 

m^nts.  *  ’  ’  The  amount  of  backscattering  would  Increase  with  both  the 

roughness  of  the  patches  and  number  of  patches. 

A  number  of  authors  have  reported  that  backscattering  levels  from  the 
undersurface  of  sea  ice  are  anomalous,  that  is,  not  the  scattering  function 
expected  from  an  isotropic  uniformly  rough  surface.  Mellen12  noted  that 
there  are  discrepancies  between  scattering  theory  and  underice  data.  The 
theory  does  not  predict  the  observed  backscattering  functions,  and  the 
underice  spatial  roughness  spectrum  deduced  from  backscattering  functions 
is  significantly  higher  than  the  spectrum  derived  from  upward -looking, 
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narrow-beam  sonar  profiles  taken  by  submarines.  Greene13  found  at  a  pack- 
ice  site  that  backscattering  strengths  increase  as  the  angle  of  incidence 
increases.  The  scattering  from  sea  ice  is  large,  as  much  as  40  db  greater 
than  from  the  wind-blown,  ice-free  sea  surface.11  As  expected,  backscat¬ 
tering  from  the  underice  surface  tends  to  increase  with  increasing  surface 
roughness  of  a  site.1^ 

Our  data  show  two  distinct  types  of  underice  surfaces:  the  underice 
ridges  and  the  relatively  smooth  areas  between  the  ridges.  The  differences 
in  scattering  at  these  interfaces  may  be  due  to  differences  in  slope  of 
the  interface,  in  reflection  coefficient,  or  in  roughness.  Theory  shows 
the  following:  (1)  scattering  cross  sections  generally  decrease  as  the 
incident  angle  9  tends  to  90°.  The  shape  of  the  function  is  sensitive  to 
the  spectrum  of  the  roughness.  (2)  The  scattering  is  proportional  to  the 
reflection  coefficient  squared.  (3)  The  cross  section  depends  upon  the 
iftns  roughness  relative  to  the  acoustic  wavelength  and  the  cross  section  is 
asymptotic  to  a  limiting  value  for  very  short  wavelengths. 

We  now  apply  the  basic  theory  and  a  model  of  the  underside  of  rhe  ice 
(Fig.  8)  to  scattering  measurements.  At  the  smooth  surface  9  is  larger 
than  9  at  the  ridges.  The  reflection  coefficient  is  probably  smaller  at 
the  smaller  angle  of  incidence  (Langlebeu  and  Pounder 15report  that  the 
reflection  coefficient  is  0.1  rear  0°  and  increases  rapidly  for  9  greater 
than  30°).  Since  our  data  show  the  scattering  from  the  ridges  to  be  large, 
the  increase  of  roughness  and  decrease  of  9  compensates  for  any  decrease 
in  the  reflection  coefficient.  As  viewed  from  a  transducer  located  below 

a  smooth  area  (between  ice  ridges),  the  sound  beam  intercepts  more  ridges 

\ 

at  smaller  grazing  angles  and  one  would  expect  the  average  scattering  to 
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Figure  8.  Sonar  Geometry  below  sea  ice 
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increase  at  small  grazing  angles.  This  could  account  for  the  anomalous 

I  O 

scattering  function  reported  by  Greene. 

In  attempting  the  inverse  problem  of  determining  the  properties  of  the 
rough  ice  undersurface  from  backscattering  measurements  we  see  obvious 
problems.  Most  of  the  scattering  comes  from  localized  patches  associated 
with  the  undersides  of  ice  ridges.  The  average  plane  of  these  scattering 
patches  may  be  inclined  30°  or  more  relative  to  the  horizontal.  This 
information  would  have  to  be  included  in  calculating  scattering  coefficients 
&s  a  function  of  grazing  angle.  If  the  area  of  the  scattering  patch  is 
less  than  the  resolution  of  the  transducer,  then  this  must  also  be  included 
in  the  data  reduction.  An  average  reflection  coefficient  should  be  deter¬ 
mined  for  the  patch.  The  spatial  spectrum  associated  with  the  scattering 
cross  section  is  then  that  for  the  rough  patch.  Presumably  the  smooth 
areas  would  have  their  own  scattering  functions  and  spatial  roughness. 

These  difficulties  could  account  fot  the  lack  of  agreement  between  theory 
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and  data  noted  by  Mellen. 
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